Spectrum of Global Magnetorotational Instability in a Narrow 

Transition Layer 

Jesse Pino 
and 

S. M. Mahajan 

Institute for Fusion Studies, The University of Texas at Austin, Austin, Texas 78712 

pinoOmail . utexas . edu 

Received ; accepted 



-2- 



ABSTRACT 



The Global Magnetorotational Instability (MRI) is investigated for a configu- 
ration in which the rotation frequency changes only in a narrow transition region. 
If the vertical wavelength of the unstable mode is of the same order or smaller 
than the width of this region, the growth rates can differ significantly from those 
given by a local analysis. In addition, the non-axisymmetric spectrum admits 
overstable modes with a non-trivial dependence on azimuthal wavelength, a fea- 
ture missed by the local theory. In the limit of vanishing transition region width, 
the Rayleigh-centrifugal instability is recovered in the axisymmetric case, and 
the Kelvin-Helmholtz instability in the non-axisymmetric case. 



1. Introduction 



The Magn e torota tional Instability (MRI), first discussed by IVelikhovl (119591 ) and 
Chandrasekharl (119601 ). became the subject of much study when it was proposed as an 
angular momentum tra nsport mechanism in hydrodynamically stable astrophy sical systems 



by Balbus and Hawley (IBalbus fc Hawley 



1991 



1998 



Hawley fc Balbus 



1999|). Much of 



the research on MRI was done using the local approximation, which uses the implicit 
assumption that the radial variation of the equilibrium flow is sufficiently small, allowing 
for modes to be sinusoidal in both the radial and a xial directions. However, t his assumption 



can lead to misleading r esults about the existence (IMahajan fc Krishan 



( jPino fc Mahajan 



20081 ) and spectrum 



20081 ) of MRI modes. When the rotation profile changes significantly. 



globally unstable modes c an be localized to the region of greatest shear. A recent study by 



Mikhailovskii et al. 



(I2OO8I ) claimed to find 'non-local' MRI modes in a configuration with 



a step-like velocity profile. However, this mode is actually a Rayleigh type surface mode. 
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possessing the greatest growth rate in the hydrodynamic hmit. We show that a similar 
system with a narrow but finite region of velocity shear can support global MRI modes if the 
shear is stable to the Rayleigh centrifugal mode when the vertical magnetic field vanishes. 
These global modes are spatially confined to the transition region, and their properties differ 
from modes found using a local analysis. In particular, there exists a discrete spectrum of 
global modes, with growth rates that differ significantly from the local prediction if the 
vertical wavelength is comparable to the transition region width. Overstability is possible 
for modes having an azimuthal wavelength component. Our investigations examine the full 
spectrum of such modes in both the axisymmetric and nonaxisymmetric case. 



Base equations 



We consider a generic rotating equilibrium in cylindrical geometry with fiow V = rQ{r)6 
and mass density p, permeated by a uniform axial magnetic field B = B^z. We include 
gravitational and pressure terms in the ideal MHD equations, so that the sy stem remains 



relev a nt to the glo 



2002 



le.g. 



Hawlev 



Dal MRI in astrophysical settings such as accretion disks (IFrank et al. 



Wang et al. 



20031) as well as to other rota ting systems such as lab oratory experiments 



2008 



Goodman fc Ji 



I2OO2I ) and stellar core collapse (lAkiyama et al.ll2003l ). 
Since the equilibrium is uniform in the axial and azimuthal directions, we may take Fourier 
transforms in those directions (we neglect any vertical stratification or gravitational terms 
in what follows). The resulting equations for th e normal modes of Lagrangian perturbation s 



1979 



^(^)^i{k,z+m9-^t)-^ to this 



Lynden-Bell fc Ostriker 



equi librium are (IFrieman fc Rotenberg 



1960 



19671 ): 



Chanmugam 



cu^p$, - 2ipuj(V ■ V)^ - :f(^) = , 



(1) 
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where 

T{C= V(7pV-^ + (|-V)p) + V-(p^)V$, 

- V(B • (5B) + (B ■ V)(5B + ((5B ■ V)B 

+ V-(p^(V.V)V-pV(V.V)0 . (2) 

Here, 7 is the adiabatic index, and p is the pressure, and $g is the gravitational potential. 
We have ignored self-gravity effects by neglecting perturbations S^g to this potential. The 
magnetic field perturbation is 5B = V x x B). The tensor divergence in the third line is 
taken with respect to the first coordinate, i.e. V • (AB) = Vi{AiBj). 

For incompressible perturbations (V ■ ^ = 0), the radial component of Eq. ([1]) and the 
divergence condition can be reduced to: 



Fm-rVT 
ar 



2 



dmr 



ir + (3) 



^m-:i-Kr) = + — (4) 

r ar r 



Where ipT (= " V)p + B ■ 6'B) is the total perturbed (gas + magnetic) pressure, 
( = + k'^r'^, ojm = {oj — mVt) is the Doppler-shifted frequency, = piy^m ~ ^z"^!); 
where va = B^/ y/ppo is the Alfven speed. = —{p'/p){rVL^ — dr^g) is the Brunt- Vaisala 
frequency. Further reduction results in a single second order differential equation in the 
radial component of the Lagrangian perturbation. 



^3 p 
I ± 1 



m >./ 



+ 



r + 



( 2m^r 



2^3 



pkz,r 
rpN^ 



dlnr 



S ^ m 



(5) 
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2.1. Cartesian Limit 

This paper investigates global MRI when the rotation rate changes only over a small 
transition region of width centered at r = tq. When the region is narrow {d <^ tq) and 
sufficiently far from the origin, it is expected that a Cartesian analysis should suffice to 
capture the essence of the mode. Thus we drop any term which decays as 1/r or faster, 
except in the case of the azimuthal mode number (m/r is replaced by the continuous 
variable ky). Choosing coordinates x = r — tq, we take constant equilibrium rotation rates 
of = fii for X < —d, and Q = Q2 for x > d. In what follows, we take a linear change in 
the rotation rate, i.e. 

il(x] = = H -, 

^ ' 2 2d 2d' 

although the analysis does not qualitatively change if different velocity profiles are 
considered. In the frame rotating with the average angular velocity Q, the Cartesian limit 
of Eq. dS]) is: 



- (f^^) = \F^k' + kl [aA^p - - klp\r^') + k'p'g 



ir. (6) 



dx 

Here, A;^ = fc^ + kl, A = —1/2 dQ/d\nr = —tqAVL/ {2d) , and the Doppler shifted frequency 
in the transition region |x| < d is Um = OJ + 2AkyX. The total change in velocity across the 
transition region is AV = roAQ = —AAd. Exterior to the region, the rotation frequency 
and density are constant, and we have : 



C =\k'- kl I e., (7) 



4:pl/fujl2 

-^1,2 



with uji^2 = uj ± 2Akyd. This admits decaying solutions of the form : 

^1,2 = exp(-fi:i_2|x|) 
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We must take the real part of the square root to be positive so as to ensure that the modes 



are bo unded. The cartesian approximation ta . 



model (INarayan et al 



1987 



Balbus &: Hawley 



<.en h ere is similar to the shearing sheet 



19981 ). however, our model has well defined 



boundaries which admit spatially decaying modes exterior to the transition region. 



3. Step-like transition 



In this section, we show that the MRI is not present when 



step-like transition, (i.e. the limit d ^ 0). iLi fc NarayanI (120041 ) considered a similar 



Eq. 



is solved for a 



situation in the context of QPO's generated by Kelvin- Helmholtz instability at an accretion 



Mikhailovskii et al 



(12008 



disc-magnetosphere boundary layer. This limit was also taken by 
hereafter M08), but resulted in the incorrect labeling of the instability as MRI. Since the total 
pressure must be the same on both sides of the interface, the c omponent of the Lagrang ian 



pertu rbation normal to the interface, must be continuous (iLynden-Bell &: Ostrikei 



19671 ). but C,r can have a jump. Integrating across the narrow transition region, we arrive at: 



F2K2 + FiKi = Ap{k'g - kiroQ') - Akipron{An) 



(9) 



For axisymmetric modes, k = and cui 2 — ^ uj- Since t he frequency then app ears only 



Chandrasekhar 



1960l ): therefore 



through u^, one can show that u"^ must be real (see, e.g. 
the frequency is purely real or purely imaginary. For an unstable root, we can take 
uj Z7, Fi^2 ^(7^ + k'lv'ji^. For the remainder of the present paper, we will take the 
density constant across the boundary. Then Eq. ([9]) simplifies to: 



(10) 



This result is equivalent to the mode equation (26) of M08, in the limit of small rotation 
shear. We first note that instability is only possible for Q2 < ^i- Solving for 72, we find an 
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unstable branch: 

Taking the hmit — ^ 0, we arrive at the hydrodynamical result: 

7' = 2^2 (^1^1 + kiri{n^ - ^]2)V^' - 1) 

Near marginal stability, with 7 <^ fi, we find 7 = kzrQ{VLi — VL2)/2 = —kz{AV)/2. 
Furthermore, since 7^ in Eq. (ITT!) decreases with increasing va, the growth rate for a given 
wavelength is maximum when the field is absent. In the hydrodynamic limit, the instability 
exists at all wavelengths, the growth rate increases without bound. The addition of a 
vertical magnetic field acts as a surface tension, reducing the growth rate of all wavelengths, 
and stabilizing modes with wavenumbers above kz^cru = 2Q\AV\/v\. 

In the fully azimuthal case kz = 0, n = ky = k, and Eq. simplifies to: 

(uj + ky{AV)/2)^ + iuj- ky{AV)/2)y = (12) 
This equation admits an imaginary solution with gr owth rate 7 = ia; = kJ AV\/2. This is 



Chandrasekhar 



1961 



jlOl). We note 



the standard Kelvin-Helmholtz (KH) instability (c.f. 
that there is no centrifugal effect; the instability is due solely to the abrupt velocity change 
across the boundary. The growth rate of this instability is the same as for the axisymmetric 
Rayleigh instability if kz is replaced by ky, however, the vertical magnetic field has no 
stabilizing effect in the pure azimuthal case. When both ky and kz are nonzero, we again 
find that the magnetic field has only a stabilizing effect. If the density is not constant 
across the jump, add itional destabilization can arise due to the Rayleigh- Taylor instability 



(iLi fc Narayan 



200J) 



In contrast to these results, the ideal magnetorotational instability occurs when a 
weak magnetic field destabilizes a plasma with a negative gradient in Q, because the frozen 



8 



flux condition of MHD provides a coupling between adjacent fluid elements. For a given 
wavelength, there is a critical magnetic field strength which maximizes the growth rate of 
the MRI. This leads us to conclude that, contrary to the assertions in M08, treating a sharp 
velocity change in a rotating plasma as a step-like transition admits only magnetically 
stabilized hydrodynamical instabilities, but not the MRI. 



The step-like analysis that results in Eq. (Q fails to capture the MRI for two reasons. 
It makes the shear essentially infinite, and makes the implicit assumption that C,r is constant 
throughout the transition region. A more robust treatment of the problem is to retain the 
finite transition region of width 2d and solve for on the interior. Normalizing lengths to 
d and all frequencies to Q, Eq. ([6]) becomes: 



Where x runs from -1 to 1, K = kd is the unitless wavenumber, and A = A/Cl. In 
order to maintain consistency with the cartesian limit, we must maintain d <^ tq. If 
A = ro(i72 — Qi)/{2dri) is not too large, we can still treat the whole system as rigidly 
rotating with a small additional shear flow. Corrections to this model are 0{A'^d'^ /r^) <^ 1. 
Global modes can now be found by matching ^ and ^' at the left and right boundaries to 
the exterior solutions. 



If there is no azimuthal component to the mode wavelength {K = = kzd), the only 
position dependence is in the density. In the constant density case, the mode equation can 



4. 



Finite-Width Transition 




(13) 



4.1. Axisymmetric modes 
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be solved analytically, 



Here u;^ = fc^fA = K^B^/ (Od^/iop) is the normalized Alfven frequency. The eigenmode 
solutions are sinusoidal or exponential, depending on the sign of K^. Both external 
boundary conditions can be simultaneously satisfied only when is positive. The local 
limit is achieved when oo, and gives an upper bound to the most unstable mode. 

Provided A < 1, the maximum local growth rate is 7 = A, occuri ng at an Alfven frequeii cy 



of ujyi = a/A(2 — A). These are characteristics of the local MRI ( Balbus fc Hawley 



19981). 



The instability is cut off for fields above = 2v A . We also see that for small magnetic 



field strength, the growth rate scales as 7 = lua^/ A/{1 — A). If 1 < A < 2, the system is 
hydrodynamically unstable, but small magnetic field strengths will amplify the instability. 
Above A = 2, the instability becomes the magnetically stabilized Rayleigh-centrifugal type, 
as in the previous section. 

For finite K^, we can approximate solutions by taking rigid wall solutions {^r = 0) at 
the boundaries. Thus Kr = mT/2 for integer n. The eigenfunctions are real, and the mode 
number label denotes the number of times that solution crosses the origin. The system 



exhibits Sturmian behavior, that is, higher values of n are more stable (IGoedbloed fc Poedts 



20041 ). There is an unstable mode for every n up to a critical Umax which satisfies the 



marginal stability equation: 

When Kz is large, there are many discrete modes which satisfy the global dispersion 
relation. The number of modes increases without bound as increases with fixed loa- 
However, for fixed magnetic field strength va^ the Alfven frequency will also increase with 
reducing rimax, as well as the growth rate of the fastest growing mode. The critical 
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parameter dictating the importance of these global effects is the ratio of the Alfven transit 
time across the gap to the rotation rate, since for the fastest growing mode, Qd/vA ~ K^. 
Even if the field is weak compared to the local rotation speed va ^ '"o^j global effects can 
be important in the d <^ ro limit. The number of unstable modes is reduced for small Kz. 
Additionally, we have assumed infinite boundaries in the z direction, but in real systems, 
the vertical system size imposes a minimum = 2Tid/ L^. 

A shooting and matching code was used to find the discrete solutions that satisfy the 
boundary conditions that match Eq. ([H]). Figure [T] shows the most unstable global modes 
as a function of uja for fixed Kz = 0.1, 0.5, 1, 5, 10. Note that the growth rate vanishes for 
vanishing magnetic field, one of the primary traits of the MRI. For smaller values of kzd, 
the maximum growth rate and the cutoff frequency are much lower than in the local limit. 
Figure [2] shows the growth rate of all unstable axisymmetric modes for when Qd/vA = 3.0. 
As Kz (equivalently, ua) is increased, the modes with smaller growth rate begin to be cut 
off, until there is only one mode remaining. 



The introduction of the azimuthal wavenumber removes the Hermiticity of the problem, 
and the frequency is no longer guaranteed to be purely real or imaginary, i.e., overstability 
is possible. However, if we take K to be large with respect to the radial change in ^, and 
restrict our analysis to small x, we can neglect this and obtain the local dispersion relation: 



4.2. Nonaxisymmetric modes 



4-2.1. Local theory 



(to^ - Lulf + {4:A{uj^ - ul) - Atu^) cos^ 6 



(16) 



where cos 6' = Kz/K. This yields a growth rate of 



7^ = -ul - 2{Q - A) cos^ 9±2n cos 9 Jul + (^] - A)2 cos2 9, 
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Most unstable axisymmetric mode, A=.75 




Fig. 1. — The most unstable mode (dotted lines) as a function of lva — k^VA for A — 
0.75, Kz = kzd = 0.1,0.5,1,5,10. When is small, the maximum growth rate and the 
cutoff frequency are much lower than in the local limit (solid line). 
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Fig. 2. — Growth rates of all 5 global modes for Vtd/vA = 3.0, A = 0.75. In this case 
Kz — 3.0 uja/^- Global effects are more important for the smaller Alfven frequencies. 
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which has a maximum of 7 = A cos 6, at uja = a/A(2 — A) cosO. The instabihty is cut off 
above 00 a = 2^/Acos6. The only prediction of the local theory is that the addition of an 
azimuthal component of the wavelength is stabilizing compared to the axisymmetric case. 



4-2.2. Global analysis 

The local approximation is a poor one in the non-axisymmetric all dependence 

on position is neglected. This approach fails to capture overstability, and sheds no light 
on the effect of finite ky on the spectrum of MRI. In particular, there is the introduction 
of a continuum of Alfven singularities on the real frequency axis (see, for example 



Hasegawa &: Uberoi 



1982 



for a discussion of the mathematical structure of this continuum). 
For every value of —1 < x < 1, there are two values of u which satisfy {uj + 2AKyxy = 
Thus there are two singularity regions, oj G {icj^ — lAKy.^ojA + 2AKy}. When 
uja < 2AKy, the two regions overlap. 

We have investigated the Ky 7^ mode equation (fT3|) with a sh ooting and m atching 



19971 ). Figure 



code. All complex roots were found using a Lemur-Schur algorithm (lActonl 
[3] shows the dependence on the magnetic field strength of the most unstable axisymmetric 
modes for fixed Ky and Kz, but at a variety of different pitch angles tan^ = Ky/K^. We 
have taken a small value of kd = 0.5. We first notice that the cutoff Alfven frequency 
is the same for different angles, contrary to the predictions of the local analysis. When 
the magnetic field is low, the instability is suppressed compared to the axisymmetric case. 
There is also a real part of the frequency, which indicates that the mode is co-rotating with 
the plasma at some point x 0. As the magnetic field is decreased, this corotation point 
approaches the boundary, where Ur = ±2AKy = ±2AK sm9. 

More can be learned by examining the structure of all unstable modes for increasing 
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Ky. Figure H] shows a representative example with fixed = 5.0 and uja = 1-0. In 
the axisymmetric hmit, there are 5 purely growing modes for these parameters. As Ky 
increases, the mode frequencies remain purely imaginary, although the eigenf unctions 
become complex (we retain the mode number ordering, although n no longer refers to the 
number of nodes). The growth rates of the some of the modes decrease, while those of 
others increase. At Ky ~ 0.1, the n = and n = 1 modes have the same growth rate, and 
above that, both modes have the same growth rate, but they now have a real part. Due to 
the symmetry of the mode equation, if ujr + "i^^i are solutions, then so are ±u!r ± iuJi. Here 
we only plot unstable roots with positive real part (dashed lines in Fig. H]). As Ky increases 
further, the real part increases linearly and the growth rate decreases. The growth rate of 
the most unstable global mode agrees fairly well with the local theory in the axisymmetric 
case, but the discrepancy grows rapidly as Ky increases. 

At larger Ky, higher order mode pairs merge as well. At Ky ~ 0.28, a new purely 
growing mode (n = 5) becomes accessible, and eventually others are created as well. When 
Ky ~ 1.0, the n = 6 complex mode merges with the real Alfven continuum. This trend 
continues until all instabilities are cut off above K,, ~ 4.3. 



4.3. Effective Potential Formulation 

The behavior of the nonaxisymmetric modes for varying Ky can be understood in 
the following manner. If we make the transformation y = (u;^ — uj]^^r then Eq. (fT3ll 
becomes: 

/ , + 4^2 2 _ 4^2^(^2 _ ^2 ) N 

y'\x) =[K^ " yi^) = ^)?/(^)' (17) 

recalling that uJm = uj + 2AKyX. For a given set of parameters, there is an effective 
(complex) potential for the perturbation. If the y associated with a given potential can 
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Fastest Growing Mode, K=0.5 



7: 

0.5 r 




Fig. 3. — (Color online) Growth rates (solid) and real components (dashed) of the frequency 
for several angles (tan^ = KyjK^^ with K = 0.5, A = 0.75. For smaller ua, the real part 
approaches the Alfven singularity value at the boundary, uj = 2^47^ sin 6*, and the growth 
rate vanishes. For larger ua, there are still pure imaginary modes which have the normal 
MRI dependence, but these decrease in growth rate as the angle is increased. 
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satisfy both boundary conditions, then that perturbation is an eigenmode of the system. 
The real part of the potential must therefore be sufficiently negative over enough of the 
transition region if the \y\ is to decay on both sides of the region (since |y|V|y| must change 
sign). Axisymmetric MRI eigenmodes have a V which is a negative real constant over the 
region. Decreasing 7 = '^{00) will make V more negative, so smaller growth rates correspond 
to more spatial oscillations, in agreement with the previous section. In the nonaxisymmetric 
case, if u is purely imaginary, ^{V) is symmetric about a; = (5^(1^) is antisymmetric). As 
long as 7 7^ 0, there is no Alfven singularity, but the potential changes the fastest near the 
points of closest approach, x = (— ±uJA)/{'2AKy). The real part of V resembles a positive 
resonance curve near these points, with the growth rate 7 directly related to the width of 
the resonance. These 'near-resonances' act as barriers, and the modes are evanescent inside 
them. When Ky is very small, these points lie far outside the transition region, and a given 
eigenmode remains similar to its axisymmetric counterpart, a seen in Figure [5^. As Ky 
increases, the resonances move closer to the boundary. Each mode at this point is related 
to either the left or right resonance, so it becomes necessary for the potential to shift, 
so as to keep only one resonance in the transition region (Figure Et). For an overstable 
mode, the ^(V) is symmetric around Xq = —u!r/{2AKy). As the growth rate decreases and 
Ur increases, the mode moves further over and becomes more oscillatory, approaching an 
Alfven continuum mode at the boundary x = ±1. 



5. Conclusions 

Global studies of the magnetorotational instability tend to depend strongly on specific 
boundary conditions. When the rotation rate changes only in a small area, however, velocity 
shear can give rise to an effective potential which spatially confines global eigenmodes and 
reduces dependence on the boundary conditions. We have examined an idealized case where 
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y, 0) 




05 1.0 1.5 



Fig. 4. — (Color online) All unstable global modes for = 5.0, A = 0.75, u;^ = 1.0. As 
Ky is increased, pairs of purely growing modes become closer together in growth rate. At 
the point of merger, the frequencies have a real part (dashed lines) which starts at zero and 
increases linearly as Ky increases. Near Ky = 0.28, a new purely growing mode is created. 
At a critical Ky for each mode, the growth rate vanishes and the mode merges with the 
(stable) Alfven continuum. The dotted line gives the maximum growth rate predicted by 
the local theory. Although the n = mode is close to this limit in the axisymmetric case, 
the global growth rates are significantly smaller for Ky ^ 0. 
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0.52722 i Effective Potential 




Fig. 5. — (Color online) a) Eigenmode structure of n = 2 mode {uo = 0.5272i) for = 5, 
Ky = 0.2, uja/^ = 1., and A = 0.75. The solid (dashed) line is the real (imaginary) part 
of C,r)- b) Effective potential for the mode shown in (a). Since Ky is small, V is changed 
little from the axisymmetric case of a negative real constant, c) For Ky = 0.75, the same 
mode now has a complex frequency {u = 0.4956 + 0.2245i), and has thus shifted its center 
to xq = —uJr/{2Aky) = —0.472. d) The effective potential of the mode shown in (c) displays 
the near- resonance barrier on the right. Since the mode has shifted to the left, the other 
peak is outside the transition region. 
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Q{r) changes sharply between two area s of rigid rot a tion; such a configuration has potential 



application to both astrophysical (e.g. 



Heger et al.l (120001 ) . where sharp rotation gradients 



arise in models of stellar core collapse) as well laboratory settings (e.g. Ekman fiow). When 
the velocity change is treated as a step-like transition as done in M08, the axisymmetric 
instability is most unstable in the limit of vanishing magnetic field- we thus conclude that 
this is a Rayleigh type surface mode. The step-like transition is Kelvin-Helmholtz unstable 
as well in the non-axisymmetric regime. When the shear region has non-vanishing thickness, 
the unstable Rayleigh mode persists when A = —l/2dlnQdlnr > 1. However, when the 
shear is moderate, {A < 1) the fiow is hydrodynamically stable, but is destabilized by the 
addition of a small vertical field. This is the hallmark of the MRI. There is a spectrum of 
purely growing global axisymmetric instabilities. If the vertical wavelength is comparable 
to the transition region width, the growth rates of the most unstable global modes differ 
significantly from the predictions of the local theory. 

That the MRI results from t he de stabilization of the slow Alfven wave by rotation 



was shown in iBalbus fc Hawleyi (119981 ). In the present analysis, we have shown that 
there is a spectrum of non-axisymmetric unstable and overstable global modes which 
connect smoothly to the Alfven continuum. For a given axial wavenumber kz, there is a 
complex dependence on the azimuthal wavenumber ky. There exists a minimum azimuthal 
wavelength below which perturbations are stable. The interaction of these modes can have 
important consequences for turbulence in the nonlinear regime-this will be the subject of 
future work. 

The authors thank Dr. Richard Hazeltine and Dr. J. Craig Wheeler for useful 
discussions. 
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